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Abstract 

The performance in finite time of a discrete heat engine with internal friction 
is analyzed. The working fluid of the engine is composed of an ensemble 
of noninteracting two level systems. External work is applied by changing 
the external field and thus the internal energy levels. The friction induces 
a minimal cycle time. The power output of the engine is optimized with 
respect to time allocation between the contact time with the hot and cold 
baths as well as the adiabats. The engine's performance is also optimized 
with respect to the external fields. By reversing the cycle of operation a 
heat pump is constructed. The performance of the engine as a heat pump is 
also optimized. By varying the time allocation between the adiabats and the 
contact time with the reservoir a universal behavior can be identified. The 
optimal performance of the engine when the cold bath is approaching absolute 
zero is studied. It is found that the optimal cooling rate converges linearly to 
zero when the temperature approaches absolute zero. 

I. INTRODUCTION 

Analysis of heat engines has been a major source of thermodynamic insight. The second 
law of thermodynamics resulted from Carnot's study of the reversible heat engine [0. Study 
of the endo-reversible Newtonian engine |§ began the field of finite time thermodynamics 
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]. Analysis of a virtual heat engine by Szilard led to the connection between thermo- 
dynamics and information theory [7||§. Recently this connection has been extended to the 
regime of quantum computation ||. 

Quantum models of heat engines show a remarkable similarity to engines obeying macro- 
scopic dynamics. The Carnot efficiency is a well established limit for the efficiency of lasers 
as well as other quantum engines [JlQHlljl . Moreover, even the irreversible operation of quan- 
tum engines with finite power output has many similarities to macroscopic endo-reversible 
engines 
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It is this line of thought that serves as a motivation for a detailed analysis of a discrete 
four stroke quantum engine. In a previous study p0 |, the same model served to find the 



limits of the finite time performance of such an engine but with the emphasis on power 
optimization. In that study the working medium was composed of discrete level systems 
with the dynamics governed by a master equation. The purpose was to gain insight into 
the optimal engine's performance with respect to time allocation when external parameters 
such as: the applied fields, the bath temperatures and the relaxation rates were fixed. 

The present analysis emphasizes the reverse operation of the heat engine as a heat pump. 
For an adequate description of this mode of operation inner friction has to be a consideration. 
Without it the model is deficient with respect to optimizing the cooling power. Another 
addition is the optimization of the external fields. This is a common practice when cold 
temperatures are approached. With the addition of these two attributes, the four stroke 
quantum model is analyzed both as a heat engine and as a refrigerator. 

Inner friction is found to have a profound influence on performance of the refrigerator. 
A direct consequence of the friction is a lower bound on the cycle time. This lower bound 



excludes the non-realistic global optimization solutions found for frictionless cases ]2D] where 
the cooling power can be optimized beyond bounds. This observation, has led to the sug- 
gestion of replacing the optimization of the cooling power by the optimization of the cooling 
efficiency per unit time Including friction is therefore essential for more realistic 

models of heat engines and refrigerators with the natural optimization goal becomes either 



the power output or the cooling power. The source of friction is not considered explicitly 
in the present model. Physically friction is the result of non-adiabatic phenomena which 
are the result of the rapid change in the energy level structure of the system. For example 
friction can be caused by the missalignement of the external fields with the internal polariza- 
tion of the working medium. For a more explicit description of the friction the interactions 
between the individual particles composing the working fluid have to be considered. The 



present model is a microscopic analogue of the Ericsson refrigeration cycle |25| where the 
working fluid consists of magnetic salts. The advantage of the microscopic model is that 



the use of the phenomenological heat transfer laws can be avoided |T6[ . The results of the 



present model are compared to a recent analysis of macroscopic chillers ||27|| . In that study, 
a universal modeling was demonstrated. It is found that the discrete quantum version of 
heat pumps has behavior similar to that of macroscopic chillers. 

There is a growing interest in the topic of cooling atoms and molecules to temperatures 
very close to absolute zero [p%[| . Most of the analysis of the cooling schemes employed are 
based on quantum dynamical models. New insight can be gained by employing a thermo- 
dynamic perspective. In particular the temperatures achieved are so low that the third 
law of thermodynamics has to be considered. The discrete level heat pump can serve as 
a model to study the third law limitations. The finite time perspective of the third law is 
a statement on the asymptotic rate of cooling as the absolute temperature is approached. 
These restrictions are imposed on the optimal cooling rate. The behavior of the optimal 
cooling rate as the absolute temperature is approached is a third law upper bound on the 
cooling rate. The main finding of this paper is that the optimal cooling rate converges to 
zero linearly with temperature, and the entropy production reaches a constant when the 
cold bath temperature approaches absolute zero. 
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II. BASIC ASSUMPTIONS AND FORMAL BACKGROUND FOR THE HEAT 

ENGINE AND THE HEAT PUMP 



Heat engines and heat pumps are characterized by three attributes: the working medium, 
the cycle of operation, and the dynamics which govern the cycle. Heat baths by definition 
are large enough so that their temperatures is constant during the cycle of operation. The 
heat engine and the heat pump are constructed from the same components and differ only 
by their cycle of operation. 



A. The Working Medium 

The working medium consists of an ideal ensemble of many non-interacting discrete level 
systems. Specifically, the analysis is carried out on two-level systems (TLS) but an ensemble 
of harmonic oscillators |[0) would lead to equivalent results. 

The TLS systems are envisioned as spin- 1/2 systems. The lack of spin-spin interactions 
enables the description of the energy exchange between the working medium and the sur- 
roundings in terms of a single TLS. The state of the system is then defined by the average 
occupation probabilities P + and P_ corresponding to the energies \oj and where to is 

the energy gap between the two levels. The average energy per spin is given by 

E=p + .fa)+p-(-y) (2.D 

The polarization, S, is defined by 

S = \{P + ~ P-) , (2-2) 

and thus the energy can be written as E = ujS. Energy change of the working medium can 
occur either by population transfer from one level to the other (changing S) or by changing 
the energy gap between the two levels (changing uj). Hence 

dE = Sdu + oodS . (2.3) 
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Population transfer is the microscopic realization of heat exchange. The energy change due 
to external field variation is associated with work. Eq. (|2.3| ) is therefore the first law of 
thermodynamics: 

DW = Sdu ■ DQ = oodS . (2.4) 
Finally, for TLS the internal temperature, T', is always defined via the relation 

Note that the polarization S is negative as long as the temperature is positive. 



B. The Cycle of Operation 

1. Heat engines cycle 

The cycle of operation is analyzed in terms of the polarization and frequency (S, u). A 
schematic display is shown in Fig.(^J) for a constant total cycle time, r. The present engine 
is an irreversible four stroke engine ]20|] resembling the Stirling cycle, with the addition of 
internal friction. The direction of motion along the cycle is chosen such that net positive 
work is produced. 
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FIG. 1. The heat engine with friction in the u ,S plane. Th is the hot bath temperature, Th is 
the time allocation when in contact with the hot bath. T c and r c represent the temperature and 
time allocation for the cold bath. r a represents the time allocation for compression (field change 
from ujf) to uj a ) and for expansion. The area A,B, C 1 , D 1 is the positive work done by the 
system, while the areas C, C 1 , Si, S3, and D 1 , D, S4, S2 represent the negative work done by the 
system. 

The four branches of the engine will be now briefly described. 

On the first branch, A — > B, the working medium is coupled to the hot bath of temper- 
ature Th for period Th, while the energy gap is kept fixed at the value Ub- The conditions 
are such that the internal temperature of the medium is lower than Th. In this branch, 
the polarization is changing from the initial polarization 52 to the polarization Si. The 
inequality to be fulfilled is therefore: 



Since u is kept fixed, no work is done and the only energy transfer is the heat Ub(Si — S 2 ) 
absorbed by the working medium. 

In the second branch, B —>■ C the working medium is decoupled from the hot bath for a 
period r tt , and the energy gap is varied linearly in time, from cj& to u a . In this branch work is 
done to overcome the inner friction which develops heat, causing the polarization to increase 
from Si to S3 (Cf. Fig. [I]). The change of the internal temperature is the result of two 
opposite contributes. First lowering the energy gap leads to a lower inner temperature for 




(2.6) 



6 



constant polarization S. Second increase in polarization due to friction, leads to an increase 
of the inner temperature for fixed u. The inner temperature T' at point C might therefore 
be lower or higher than the initial temperature at point B. 

The third branch C — > D, is similar to the first. The working medium is now coupled to 
a cold bath at temperature T c for time r c . The polarization changes on this branch from S3 
to the polarization S4. For the cycle to close, £4 should be lower than £2. At the end of the 
cycle the internal temperature of the working medium should be higher than the cold bath 
temperature, T" > T c , leading to: 

^4 > -^tanhf-^-) . (2.7) 



2 \2k B T c 

Since S4 < Si (Fig. §, it follows from Eq. ( |2j) and Eq.( g/7|), that: 

(i) > m ™ 

Inequality ( |2.8| ) is equivalent to the Carnot efficiency bound, from Eq. (|2.8|) one gets: 

1 " (S) < 1 " (I) = (2 ' 9) 

The present model is a quantum analogue of the Stirling engine which also has Carnot 's 
efficiency as an upper bound. 

The polarization S changes uni-directionally along the 'adiabats' due to the increase of 
the excited level population as a result of the heat developed in the working fluid when work 
is done against friction, irrespective of the direction of the field change. 

The fourth branch D — > A, closes the cycle and is similar to the second. The working 
medium is decoupled from the cold bath. In a period r& the energy gap is changing back to 
its original value, u>t,. The polarization increases from S4 to the original value S^- 
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TABLES 



TABLE I. Work and heat exchange along the branches of the heat engine with friction 



branch 


work+[work against friction] 


heat 


A -> B 





Ub(Si - S 2 ) 


B^C 


(u a - w b )(5i + a 2 /(2r a )) + [ a 2 {to a + w 6 )/(2r«) ] 





C -^D 





^ a ((5 2 - Si) - a 2 (l/r a + 1/Tfc)) 


D -> A 


(u; fe -w a )(5 2 -a 2 /(2r fe )) + [ a 2 (w a + w 6 )/(2r 6 ) ] 






Refrigerator cycle 

The purpose of a heat pump is to remove heat from the cold reservoir by employing 
external work. The cycles of operation in the (S, u) plane is schematically shown in Fig. |2|, 




CO 



FIG. 2. The cycle of operation of the heat pump. Left: without friction. S'f 9 is the hot bath 
equilibrium polarization. S^ 9 is the cold bath equilibrium polarization. The area enclosed by 
D,C,S2 Si is the heat absorbed form the cold bath. The area enclosed by DCBA is the work done 
on the system. Right: with friction. The area enclosed by D,C,S2,Si is the heat absorbed form the 
cold bath. The work on the system is the area defined by the rectangles B B 1 S 2 S3 and B 1 A 1 D 
C and A 1 A S 4 S v 
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The cycle of operation resembles the Ericsson refrigeration cycle ||25|| . The differences 
are in the dynamics of the microscopic working fluid which are described in subsection [11 C . 
The work and heat transfer for the heat pump is summarized in Table [TT]. 

The four branches for the heat pump become: 

In the first branch, D — > C, the working medium is coupled to the cold bath of tem- 
perature T c for time r c , while the energy gap is kept fixed at the value uo a . The conditions 
are such that the internal temperature of the medium is lower than T c during r c . Along 
this branch, the polarization changes from the initial polarization S\ to the polarization S2. 
Since 10 is kept fixed, no work is done and the only energy transfer is the heat u) a {S2 — Si) 
absorbed by the working medium. On this branch: 

S 2 < -Jtanhf-^-) . (2.10) 



2 \2k B T c 

In the second branch, C — > B the working medium is decoupled from the cold bath, and 
the energy gap is varied. In the frictionless case the polarization S2 is constant (Left of Fig. 
0). The only energy exchange is the work done on the system ( Table [TJ). When friction is 
added the polarization is changing from S2 to S3 in a period r a . The energy gap changes 
from uj a to Ub (Right of Fig. g), according to a linear law. In addition to work, heat is 
developing as a result of the inner friction ( Table |I|). 

The third branch B — > A, is similar to the first. The working medium is coupled to the 
hot bath at temperature T h , for time Th, keeping the energy gap ui^ fixed. In this branch the 
polarization changes from S2 to Si in the frictionless case, and from S3 to S4 when friction 
is added. The constraint is that the internal temperature of the working medium should be 
higher than the hot bath temperature during the time r h , T' > T h , leading to the inequality 
(Fig. 0), 

Si > S 4 > - - tanh ( . (2.11) 



2 \2k B T h , 

therefore S2 > Si. From Eqs. ( |2.1(J| ) and ( |2.11| ), the condition for the interrelation between 



the bath temperatures and the field values becomes: 

9 



which is just the opposite inequality of the heat engine, (Eq. |2.8| ). In the heat pump work is 
done on the working fluid and since no useful work is done Carnot's bound is not violated. 

The fourth branch A —>■ D, closes the cycle and is similar to the second. The working 
medium is decoupled from the cold bath, and the energy gap changes back, during a period 
r b to its original value, u b . 

The results are summarized in Table ||. 



TABLE II. Work and heat exchange along the branches of the heat pump without /with friction. 



branch 


frictionless work/work+[work against friction] 


heat 


D^C 





w«(5 2 - Si) 


C^B 


(u b - UJ a )S 2 

(uj b - uj a )(S 2 + cj 2 /{2T a )) + [a 2 {uj a + u b )/(2r a )} 





B^A 





uj b (Si - s 2 ) 

co b {(St - S 2 ) - a 2 (l/r a + 1/n)) 


A^D 


(u a - u b )Si 

K - oj b ) (Si - a 2 /{2r b )) + [a 2 {u a + u b )/{2 n )\ 
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C. Dynamics of the working medium 



The dynamics of the system along the heat exchange branches is represented by changes 
in the level population of the two-level-system. This is a reduced description in which 



the dynamical response of the bath is cast in kinetic terms [ 18[ . Since the dynamics has 



been described previously |2(| only a brief summary of the main points is presented here, 
emphasizing the differences in the energy exchanges on the 'adiabats'. 

1. The dynamics of the heat exchange branches 

The dynamics of the population at the two levels, P + and P_, are described via a master 
equation 

dt , (2.13) 

it = hP + - M>- 

where and fcj are the transition rates from the upper to the lower level and vice versa. 
The explicit form of these coefficients depend on the nature of the bath and the system bath 
coupling interactions. The thermodynamics partition between system and bath is consistent 



with a weak coupling assumption Temperature enters through detailed balance. The 
equation of motion for the polarization S obtained from Eq. ( |2.13| ) becomes: 

-P = -T(S - S eg ) (2.14) 
dt 

where 



h + fef (2.15) 



and 



S« = - l - k -^J?± = -Itanh (2.16) 
2k l + fc T 2 \2k B TJ v ; 

where S eq is the corresponding equilibrium polarization. It should be noticed that in a TLS 

there is a one to one correspondence between temperature and polarization thus internal 

temperature is well defined even for non-equilibrium situations. 
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The general solution of Eq (|2.14j) is 



S(t) = S eq + (5(0) - S eq )e~ vt . (2.17) 



where S(0) is the polarization at the beginning of the branch. 
From Eqs. fl2.14p and Q2.16| ) the rate of heat change becomes: 



Q = uS (2.18) 



See also [16 . 



For convenience, new time variables are defined: 

x = e~ TcTc , y = e~ VhTh (2.19) 

These expressions represent a nonlinear mapping of the time allocated to the hot and cold 
branches by the heat conductivity V. As a result, the time allocation and the heat conduc- 
tivity parameter become dependent on each other. 

Figure [IJ and ^ show that the friction induces an asymmetry between the time allocated 
to the hot and cold branches since more heat has to be dissipated on the cold branch. 



2. The dynamics on the 'adiabats' 

The external field u and its rate of change Co are control parameters of the engine. For 
simplicity it is assumed that the field changes linearly with time: 

oj{t) = Cot + u(0) (2.20) 

Rapid change in the field causes non-adiabatic behavior which to lowest order is propor- 
tional to the rate of change Co. In this context non-adiabatic is understood in its quantum 
mechanical meaning. Any realistic assumption beyond the ideal non-interacting TLS will 
lead to such non-adiabatic behavior. It is therefore assumed that the phenomena can be 
described by a friction coefficient a which forces a constant speed polarization change S: 
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e> 2 



S = ^ (2.21) 

where £' is the time allocated to the corresponding 'adiabat'. Therefore, the polarization as 
a function of time becomes: 

S(t) = 5(0) + (^)\ (2.22) 



M) 

where t > , t < t' . A modeling assumption of internally dissipative friction, similar to 
Eq. (|2.21 ), was also made by Gordon and Huleihil ( pB[). Friction does not operate on the 
heat-exchange branches, there is no nonadiabtic effect since the fields u a and u>b are constant 
in time. The irreversibilities on those branches are due to the transition rates (r) of the 
master equation. 

From Fig. ([!]), Eq. ( |2.4| ), and Eq. ( |2.22| ) the polarization, for the B — > C branch of the 
heat engine becomes: 



a 2 



S C = S 3 = S 1 + (_J . (2.23) 
The work done on this branch is: 

W BC = [ Ta DW = r Sudt = {u a - u b )[Si + (2.24) 



The heat generated on this branch in the working fluid, which is the work against the friction, 
becomes: 

Ta f Ta ri (T 2 (uJ a + UJ b ) 



Qbc = / DQ = / uSdt = v ; ' 01 (2.25) 
Jo Jo 2r a 

This work is dependent on the friction coefficient and inversely on the time allocated to the 

'adiabats'. The computation for the other branches of the heat engine and heat pump are 

similar. 



3. Explicit expressions for the polarizations imposed by the closing of the cycle. 



By forcing the cycle to close, the four corners of the cycle observed in Fig. [I] are linked. 
Applying Eq. ( [2.171) leads to the equations: 
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51 = S 2 y + S^(l - y) 
S 3 = S 1 + £ 

S A = S 3 x + S?(l-x) 

5 2 = 5 4 + ^ 

The solutions for Si, S 2 and Si — S 2 are 



^2.26) 



C _ qeq , AS e "(l-y) + a 2 yG(x) _ neq AS e "y(l-x) - a 2 yG(x) 

q _ Qeq AS eq x(l-y) + a 2 G(x) _ neq AS ee >{l~x) - a 2 G(x) 



^2.27) 



and 



5! - S 2 = (AS^)F(x, y) - a \ v)G . {x) (2.28) 

(1 - xy) 



where 



F(x,y)= (1 , AS e « = (5^ -ST) > G(x) = (x/ra + l/Ti) 

The constraint that the cycle must close leads to conditions on the polarizations S\ and 5*2 
and on the minimum cycle time T c , m in- Eqs. ( p. 27 ) shows that both Si and S 2 are bounded 



by S'h and S^ q . The minimum cycle time is obtained when the polarizations coincide with 
the hot bath polarization: Si=S2=S^ q . In this case, 7^=0, and from Eqs. (|2.19|) and Q2.28|) 
the minimum time allocation on the cold bath r c min is computed, 



(SI - - gi/n 
(S e h q - S?) + a 2 /r a 



_ w h Z^£l - u l' b C2 291 



or 



T . - i/n- W-^' I-gVl* (230) 

From this expression for r c m j n the lower bound for the overall cycle time, is obtained (The 
left of Fig. D : 

T~ ^ Train Terrain ~\~ T~ a -\- Tj> (2.31) 
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When the minimum cycle time Eq. (|2.30|) diverges, the cycle cannot be closed. This condi- 



tion imposes an upper bound on the friction coefficient a 



a < a up = J r b (S e h q - Sc q ). (2.32) 



or 

a 2 

7~b > T~b,min ^ geq g e 1^ ' (2.33) 

Closing of the cycle imposes similar constraints on the minimal cycle time under friction 
for the heat pump. The value of the polarization difference 5*2 — Si using the notation of 
Fig. |2| becomes: 

S 2 - S, = {S? - ST)F(x, y) - a2(1 ~^ )(y/T " + 1/rb) ( 2 - 34 ) 

(1 - xy) 

The minimum cycle time is calculated in the limit when r c =0, leading to S 2 —Si—S 2 q . From 
Eqs. fl2.19| ) and (|2.34 ) the minimum time allocation on the hot branch Th t mi n is computed: 

(S? - St q ) - a 2 /r b 



(S e 2 q - Sf 1 ) + 0*/T a 



[2.35) 



/r (SI - gn - a 2 In 

Th,min — L / L h l &, Q eq Q eq\ T 2/ ' 

[o 2 — oi ) -+- a /r a 

where S 2 q is point F and S{ q is point E on Fig. ||. Using r himin the lower bound for the 
overall cycle time, is computed 

T ^ Train = T~h,min ~\~ 7~ a + Tf, (2.37) 

Closing the cycle imposes a minimum cycle time for both the heat engine and the heat pump, 
which is a monotonically increasing function of the friction coefficient a. The divergence of 
T min imposes a maximum value for the friction coefficient a. 



D. Finite Time Analysis 
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1. Quantities to be Optimized. 



The primary variable to be optimized is the power of the heat engine and the heat-flow 
extracted from the cold reservoir of the heat pump. For a preset cycle time, optimization 
of the power is equivalent to optimization of the total work, while optimization of heat flow 
is equivalent to the optimization of the heat absorbed. The entropy production will also be 
analyzed. 

(1) The total work done on the environment per cycle of the Heat Engine. 

The total work of the engine, is the sum of the work on each branch: Cf. (Table | and 
Fig. 0): 



Wcyiei = j>DW = -(Wab + Wbc + Wcd + Wda) (2.38) 
which becomes: 

Wc, to i = (co b -co a )(S 1 -S 2 ) - a 2 u a (l/r a + l/r h ) (2.39) 

The negative sign is due to the convention of positive W when work is done on the system. 

Analyzing Eq. (|2.39|) , the work is partitioned into three positive and negative areas. The 
positive area (left rotation) 

W p = (oj b -oj a )(S 1 -S 2 ) (2.40) 

is defined by the points A, B, C 1 , D 1 in Fig. [II. The two negative areas (right rotation) 

W n = a 2 u a (l/r a ) + cr 2 uj a (l/T b ) (2.41) 

are defined by the points C, C 1 , Si, S3 and D 1 , D, S±, S 2 in Fig. [l|. 

The cycle which achieves the minimum cycle time r = T c ^ min , produces zero positive work 
W p = 0. The corners A and B coincide at E, and C 1 coincides with D . The negative work 
of Eq. ( |2.41| ), is defined by the corners C, D, S4, S3 and is 'cut' by the S^ q line (Cf. the right 



of Fig. The cycle has negative total work, meaning that work is done on the working 
fluid against friction. When r increases beyond T CiJn i n , S\ diverts from S 2 , becoming lower 
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than S^ q (Cf.Eq. (|2.27|) ). At a certain point, the work done against friction is exactly 
balanced by the useful work of the engine. The minimum time in which this balance is 
achieved is designated r . Its value which can be deduced from Eq. ( |2.39| ) is worked out in 
appendix [B]. 

The minimum cycle time T min is compared to To, the minimum time needed to obtain 
positive power shown in the right of Fig. ||| as a function of the friction o. Both functions 
increase with friction, but To diverges at a much lower friction parameter. Above this 
friction parameter no useful work can be obtained from the engine. The divergence of r m j n 
corresponds to a larger friction value where the cycle cannot be closed. 

—i 2.0 
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FIG. 3. Left: Minimal cycle time T m i n as a function of the inner friction parameter a for the 
heat engine. The vertical line represents the upper-bound of a. Dimensionless units are used in 
which hf, = 1 and h = 1. The parameters used are: uj a = 1794, ujf, = 4238, T c = 500, = 2500, 
r c = 1 and Th = 2. Right: Comparison between T m i n and To, the minimum cycle time for power 
production. 



When the total time allocation is sufficient, i.e. r > tq, work is done on the environment, 
and 5*i starts to increase. For long cycle times Si will approach S 1 ^ 1 , while S2 will approach 
Sl q . The constant negative area will become negligible in comparison to the positive area ( 
Fig. |). 

To study the influence of friction on the work output the polarization difference from 
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Eq. ( |2.28|) S1-S2 is inserted into the work expression Eq. ( [2.39|) , leading to 



where 



Wcyui = {u h -u a ){S^-St q )F{x,y) - W al 



2 (ub(i - y)(x/r a + i /n) u a (i-x)(i/T a + y/n)\ 

Wal = a 1 

V 1 — xy 1 — xy J 



(2.42) 



(2.43) 



W a i is the additional 'cost' due to friction and is always positive. 

The emergence of positive power V is shown in Fig. |j. For a fixed cycle time the 
optimization of work is equivalent to the optimization of power. 

The first two cycles have a cycle time shorter than tq, and therefore do not produce 
useful work. For cycle 3, r > tq and positive work is obtained when the time allocation on 
the cold bath is sufficient r c > 0.08. 
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FIG. 4. Left: Power as a function of the time allocation on the cold branch corresponding to 
the friction coefficient a = 0.005 with changing cycle times. The cycle time values are: for curve 
1, r = T m i n =0.059, for curve 2, r = 0.1 (the first two plots overlap) and for curve 3, r = 0.5. 
Other parameters are the same as in figure 3. The dashed horizontal line is the line of zero power. 
Right: The cycles corresponding to the power plots. Negative work is in blue and positive work is 
in red. Note that for cycles 1 and 2, the total area is negative and, therefore, the power output is 
negative. 



For longer total cycle times, the ratio between the negative area to the positive area 
decreases as can be seen in Fig. [5]. 



The position of the cycles in the S, lo coordinates relative to and S^ q chang es as a 
function of the cycle time. Insight to the origin of the behavior of the 'moving' cycles is 
presented in Fig. [H] of Appendix [A]. 
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FIG. 5. Left: Power as a function of the time allocation on the cold branch corresponding to 
the friction coefficient a = 0.005 with changing cycle times. The cycle time values are: for curve 
4, r = 1, for curve 5, r = 2 and for curve 6, r = 5. The dashed horizontal line is the line of zero 
power. Right: The cycles corresponding to the power plots. All the constant parameters are as in 
Fig. | 

The calculation of the total work done on the working fluid per cycle, VV°™ de3 for the 
heat pump is described in appendix (D]. See also ( Cf. Table (||) and Cf. Fig. |2[). 
(2) The heat-flow(Q F ) 

The heat-flow, Qp, extracted from the cold reservoir is: 



Qf = u a (S 2 - S x )/t 



[2AA) 



Due to the dependence of Qf only on S 2 -Si, the cycle is similar to the cycle of the heat 
engine. 

(3) The entropy production (AS U ). 

The entropy production of the universe, AS U , is concentrated on the boundaries with 
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the baths since, for a closed cycle, the entropy of the working fluid is constant. The compu- 
tational details for both the heat engine and the heat pump are shown in appendix |C[ The 
entropy production and the power have a reciprocal relation (See Fig. [T2| ). For example, 
the entropy production increases with a, while the power decreases. 
(4) Efficiency. 

The efficiency of the heat engine is the ratio of useful work to the heat extracted from 
the hot bath. 

Wcyde fricles ( CF 2 UJ a (l/T a + 1/t 6 ) \ 

Vh.e. - 7; - Vh.e. ~ 73 5^ — \ lAb ) 

^■absorbed \ ^bWl ~ ^2) J 

where Vh'e^ = (1 - u a /u b ) 

When the cycle time approaches its minimum r — > T min , the efficiency diverges: 
Vh.e. — ► — oo- The efficiency becomes positive only when r > r . Using Eq. (|2.45|) 
a bound for the efficiency is obtained: 



n ^ . frides T c ( a 2 (l/r a + 1/t 6 )\ 

The cooling efficiency of the refrigerator will be: 



^2.46) 



Qdc uJjh ~ Si) 

W W- de ^ ((^-^(^-SO + a^(l/r a + l/r 6 )) 1 ' J 



or: 



1 +1 = _J_ + 1 = ^ A | a 2 (l/r a + l/T 6 ) \ > T;/ i+( T 2 (l/r a + l/r 6 )' 



VRf COP u a V (5 2 - SO ; T c V (5 2 - 50 



leading to the expression for the efficiency: 



(2.48) 



W u , Jricles , a 2 (l/r a + l/r b ) < J> /rtdes , a 2 (l/r a + l/r b ) ( 2 ' 49 ) 

^Vh.E. H 5^ ±h VH.E. H 5^ 

For both the heat engine and the heat pump, the efficiency is explicitly dependent on time 
allocation, cycle time, and bath temperatures. 
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2. Optimization 



The performance of both the heat engine and the heat pump can be optimized with 
respect to: 

• (a) The overall time period r of the cycle, and its allocation between the hot and cold 
branches. 

• (b) The overall optimal time allocation between all branches. (This optimization is 
performed only for the heat pump.) 

• (c) The external fields, (ui a , u b ). 

(a) Optimization with respect to time allocation. 

The optimization of time allocation is carried out with the constant fields u a and u>b- 
The Lagrangian for the work output becomes: 

C(x,y,X) = W cyde + A (r + ln(x) + — ln(y) - r a - r b ) , (2.50) 

V J- c 1 h / 

where A is a Lagrange multiplier. Equating the partial derivatives of C(x,y, A) with respect 
to x and y to zero, the following condition for the optimal time allocation becomes: 

r c x ((1 - y)\S^ - S?) + a\l - y)(l/ Ta + y/n)) = 
T h y ((1 - xf{S e h q - Sf) - a 2 (l - x){x/r a + l/r b )) 

When <t = 0, the previous frictionless result is retrieved. (Optimizing the entropy pro- 
duction AS U leads to an identical time allocation to Eq. ( |2.51| )). 
Eq. ( |2.5ip can also be written in the following way: 



r c x((l-y)(l - y x max )) = T h y ((1 — x)(x max - x)) (2.52) 



where x max was defined in Eq. ( |2.29| ). The result is dependent on the time allocations of 
the 'adiabats', through the dependence of x max . 

For the special case when r c = T^, the relation between the time allocations in contact 
with the hot and cold baths becomes: 

21 



(2.53) 

For the frictionless case, this result coincides with the former frictionless onex = y, meaning 
that equal time is allocated to contact with the cold and hot reservoirs. When friction is 
added this symmetry is broken, Eq. ( |2.53| ), to compensate for the additional heat generated 
by friction the time allocated to the cold branch, becomes larger than the time on the hot 
branch. 

The Lagrangian for the heat-flow, Qp, extracted from the cold reservoir is defined in 
parallel to the Lagrangian for the total work. Substituting Th for T c , x for y and vice versa, 
also y ma x for x max , where y max was defined in Eq. ( p.35| ), one gets the optimal time allocation 
for the heat pump. 

Optimization of power with respect to time allocation as a function of the cycle time, r 
for different friction coefficients is shown in Fig. ^ (Left), together with the corresponding 
heat-flow (Middle) and the corresponding entropy production (Right). The left part shows 
that in the frictionless case the power obtains its maximum at zero cycle time with the value 
consistent with Eq. ( |2.58p . When friction is introduced, the maximum power decreases and 
is shifted to longer cycle times. The figure also shows, that for short times the work done 
by the system is negative, and as the friction coefficient a increases, the boundary between 
positive and negative power shifts to longer cycle times. In the Middle of Fig. || the heat- 
flow corresponding to the optimal power on the left is shown. The shapes of the power and 
heat flow curves are similar. The heat-flow values are always positive and larger than the 
corresponding power values. The entropy production (Right) shows that unlike the power 
curves the friction changes significantly the shape of the curves. The entropy production 
rate for the case with friction sharply decreases. The parallel graphs for the heat pump are 
similar. 
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FIG. 6. Left: The optimal power with respect to time allocation as a function of r, for different 
values of friction. Middle: The corresponding heat-flow (Qf)- The parameter values for both Left 
and Middle are: for plot 1, r a = Tf, =a = 0. For all the other plots r a = Tf, = 0.01 The a values for 
the curves from plot 2 to 6 are: 0.002, 0.005, 0.007, 0.0135, 0.02 respectively. Right: The entropy 
production rate corresponding to the optimal power on the left part of the figure. The additional 
curve is curve '1', which corresponds to a = 0, and r a = t\, = 0.01 . The parameter values for the 
other plots are : for plot 2, r a = =<j = 0, for all the other plots r a = Tf, = 0.01. The a for the 
curves from plot 3 to plot 7 are: 0.002, 0.005, 0.007 0.0135, 0.02 accordingly. 

(b) Time allocation optimization between all branches of the refrigerator 

Further optimization of the performance of the heat pump is possible by relaxing the 
assumption of constant time on the 'adiabats'. First the time allocation between the two 
'adiabats' is optimized, when r a + = 5, where 5 is a constant. Finally the time allocation 
between the 'adiabats' and the heat exchange branches, is optimized. These results are 



compared to the recent analysis of Gordon et. all. p7 



From Eqs. ( |2.44|) and (|2.34j) with constant time allocations along the heat exchange 



branches one gets for the cooling power: 

Qf =A - A^V- + t-^-) , (2.54) 

T a (0 - T a ) 

where Aq and A\ are constant functions of the parameters of the system. And on S, a double 
inequality is imposed r > 5 > the larger of [ (r — r^ m j n ); r& )m j n ], see Eq. ( |2.33| ). 
The optimal r a depends only on y and on 5. The optimal value of T a>opt becomes: 

T a ,opt = 5 ~ V , + ^ (2.55) 

i - y 
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Further optimization by changing the the value of 5, changes the cycle time r. This 
optimization step is done by numerical iteration. Typically the sum of the final optimal 
values of r a and r b is about twice their value before, and their ratio is about 0.7 of the value 
which was chosen initially. 

The next step is to study the time allocation between the 'adiabats' and the heat exchange 
branches when all other controls of the heat pump have optimal values. These controls 
include also the external fields of optimization which are described later. 

For comparison with Gordon et. all. H27| , the results of optimization are plotted in the 
1/ Qf ,1/rj plane for a fixed cycle time r. The following example demonstrates the method 
followed: First an optimal starting value for Qp was found which determines the time 
allocation control parameters, r c = 0.44221, r h = 0.31779, r a = 0.0084, r b = 0.0116 with 
a total cycle time of r = 0.78. Under such conditions Q,F,max = 2.9158 (1/ Q,F,max = 0.34296). 

Changing the time allocation between the 'adiabats' and the heat exchange branches 
changes the balance between optimal cooling power and efficiency. Denoting the sum r c +r^ 
by r ch , the ratio T h /r c by r hc , the sum r a + r b by r ab , the ratio r a /r b by r ab , time is transfered 
from r c h by small steps to r ab , while keeping the the ratios r^c and r ab constant. For each 
step the corresponding 1/Qf and I/77, are calculated as in Fig. 0. The relation between 
the reciprocal efficiency and the reciprocal cooling power shows the tradeoff between losses 
due to friction and losses due to heat transfer. Following the curve in Fig. [7|, starting from 
point A where the cooling power is optimal, resources represented by time allocation are 
transferred from the heat exchange branches to the 'adiabats', reducing the friction losses. 
At point B an optimum is reached for the efficiency. This point has been found by Gordon 
et. al. to be the universal operating choice for commercial chillers. Point B represents the 
optimal compromise between maximum efficiency and cooling power. 



24 



10 



3? 



o 














2 











12.2 



12 



11.8 



11.6 








0.1 0.2 0.3 0.4 

arbitrary time scale ,,vj< f 

FIG. 7. The relation between efficiency and cooling power for the heat pump. The parameters 
are: The constant optimal cycle time, r, =0.78; T c =51.49, T^=257.45, w a =47.699, 0^=600, r c =l, 
Th=2, o"=0.005 Left: Comparison between l/Q^(plot 1) and l/r/(plot 2) as a function of the 
allocated time transfer from the heat exchange branches to the 'adiabats'. Zero time is the optimal 
heat-flow time allocation. Right: The Universal plot for the heat pump. The starting optimal 
point in the plane of (1/Qjr, I/77), was (0.34296, 12.202), while the maximum efficiency point B 
is (0.6322,11.6379) and time allocation (T c ,T h ,T a ,T b )= (0.22721,0.16328,0.1636,0.2259). The insert 
shows the neighborhood of point A. 

Point A is located at the maximum cooling power. If more time is allocated to the heat 
exchange branches both 1/Qp and 1/rj will continue to increase as seen in the insert of Fig. 

(c) Optimization with respect to the fields. 

The values of the fields u a and ujf, are control parameters of the engine. In a spin 
system these fields are equivalent to the value of the external magnetic field applied on the 
system. They directly influence the energy spacing of the TLS. The work function W cyc i e , 
or equivalently the power (V) is optimized with respect to the fields, subject to the Carnot 
constraint: 



(2.56) 
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Optimal power is obtained by equating independently to zero the partial derivatives of 
Wcyde, or of V = Wcycie/r by varying u a and u b . In addition the optimal solutions have 
to fulfill the inequality constraints in Eq. ( [2.561) . As a result two transcendental equations 
in u a and uj b are obtained which are solved numerically. 
The two equations are: 



V-t7 } (AS- + o*/r a ) cosh 2 
fe^T ^ + ff2 /r a ) cosh 2 (^) 



(4 fc B T c ) 
1-x 



(2.57) 



(4 k B T h ) 

Where AS eq = 

^li' ^c 9 as defined in Eq. ( |2.28| ). Examining Eq. ( |2.57| ), and fixing the 
friction a, it is found that AS eq is an extensive function of order zero (intensive ) with 
respect to the quartet of variables u a ,T c ,Ub,Th. This means that scaling these parameters 
simultaneously will not change AS eq . Also x maX} and cosh 2 (^r) are ex tensive (order 
zero). The work function however, is extensive with order one (Eqs. (|2.42|) and (|2.43|) ). 
This property will be exploited in paragraph III. 

The optimization of power with respect to the fields is shown in Fig. |8] for the frictionless 
engine, as a function of the fields with fixed time allocation. A global maximum can be 
identified. 
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FIG. 8. Power for the frictionless engine as a. function of the fields {00 a ^ 

LVb) for constant bath 

temperatures, and constant time allocations. The maximum power is achieved at uj a = 1794 and 
uJb = 4239 where the bath temperatures are T c = 500, Th = 2500. 
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The heat pump optimization of Qf with respect to the fields is different and therefore 
will be presented in Section ITTI . 



The analysis for the optimization with respect to the fields for the entropy production 
AS U , is presented in appendix [C[ The optimal solution without friction(cr = 0) leads to 
AS^ in = 0. When a ^ 0, the minimum value of AS U is different from zero, and is achieved 
on the boundary of the region. 

E. Global Optimization of the Heat Engine 

Global optimization of the power means searching for the optimimum with respect to 
the control parameters cycle time, time allocation and the fields. An iterative procedure is 
used. 

The procedure is initiated by setting the optimal time allocation from the corresponding 
Lagrangian, with a = 0. The power becomes a product of two functions, one depending only 
on time the other only on the fields, and therefore, the fields can be changed independently of 
time. The optimal fields for the above time allocation are then sought. For the frictionless 
case, the overall time on the adiabats tends to zero. The optimal field values become 
independent of time. The value V = 107.501 is the short time limit in accordance with the 
equation: 

-ie<j O e <7VT T \ 1 1 /T i /T "\2 



v — (u b -u a )(s? - s?)(r c r h )/yr c + y/r h y (2.58) 

These fields are inserted into the expression with friction a ^ 0, and the new optimal 
times and fields are computed. The iteration converges after two to three steps, as indicated 



by Table III for a = 0.005. Notice that the location of the optimum is not very sensitive to 



the friction parameter. 
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TABLE III. Global optimization of power. The notation V max (uj a ,uJb) stands for fixed time 
allocations, and the notation T ,max (T a , t&) stands for fixed field values. The other parameters are: 
T c = 500, T h = 2500, r a = r 6 =0.01. T c = 1 and T h = 2 



a 


T 


V max {uJ a ,uJ b ) 


V max (r C: T h ) 






Tc/T 


0.005 


2 




84.46 


1794 


4239 


0.5999 


0.005 


1.367 




87.18 


1794 


4239 


0.5891 


0.005 


1.367 


88.68 




1719.1 


4036.31 


0.5891 


0.005 


1.347 




87.47 


1719.1 


4036.31 


0.58856 


0.005 


1.347 


88.704 




1718.16 


4033.67 


0.58856 



In Table |y], the extensive properties Eq. ( |2.57| ) are examined for k=2 and k=10 with respect 

1000, T h 



5000 for k=2 and 



to Table |T|. The temperature values will change to T c 
T c = 5000, Th = 25000 for k=10 . The results verify the analysis. 

TABLE IV. Global optimization of power. By multiplying the four values, T c , T^, uj a , Ub by k 
and searching first for optimal time allocation, then multiplying only the temperature values by k 



and searching for the optimal fields. All the notations and other parameters as in Table III 



a 


k 


r 


y {Lo a ,u}b) 


V max (T c ,T h ) 


w„ 


UJ b 


Tc/T 


0.005 


2 


1.367 




174.9 


3438.2 


8072.6 


0.58852 


0.005 


2 


1.367 


174.9 




3436.7 


8070.3 


0.58852 


0.005 


2 


1.347 




174.94 


3436.7 


8070.3 


0.58856 


0.005 


2 


1.347 


179.8 




3437.7 


8069.4 


0.58856 


0.005 


10 


1.347 


887.04 




17181.6 


40336.7 


0.58856 
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III. ASYMPTOTIC PROPERTIES OF THE HEAT PUMP WHEN THE COLD 
BATH TEMPERATURE APPROACHES ABSOLUTE ZERO. 



The goal is to obtain an asymptotic upper bound on the cooling power when the heat 
pump is operating close to absolute zero temperature. This requires optimizing the perfor- 
mance of the heat pump with respect to all control parameters. 

A. Optimization of the heat-flow Qp with respect to the fields and to the cooling 

power upper bound. 

The heat-flow, Qp extracted from the cold reservoir now becomes the subject of interest: 

Qf = u a (S 2 - S 1 )/t (3.1) 

or from Eq. fl2~3j ), 

Q F = (w /t) \{S 2 - >V)F(x,2/) _ — I (3.2) 

No global maximum for the Qf with respect to the fields is found. The derivative of Qf 
with respect to u b becomes: 

9Qf F(x,y)u a 1 



du b t Ak B T h cosh 2 



> (3.3) 



leading to the result that Qp is monotonic in u b . Under such conditions, u b is set, and the 
optimum with respect to u a is sought for. The derivative of Qf with respect to uo a becomes: 

^ = (s? - 5f ) - j^-Ay/ra + i/n) - u a 1 ua = o (3.4) 

du a (1-2/) 4 k B T c cosh 2 ^ n 

Introducing from Eq. (|3.4| ) the optimal value of (S^ — Sl q ) — (y/ T a + ^/ T b), into Eq. 
(ft.2|), leads to the optimal cooling rate: 



goptimum = F(x,y) hj 1 F(x,y) f^aV T c , 3 

T 4 k B T c cosh 2 Ak B r \T C ) cosh 2 ^fp^ 1 ' 
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Due to its extensivity, the ratio becomes a constant, while both u a and T c can 
approach zero. 

From Eg. (|3~5|) , an upper-bound for the cooling rate Qp is obtained: 



Ak B r V T £ 

From Eqs. (|3.6|) , when T c approaches zero, the cooling rate vanishes, at least linearly with 
temperature. This is a third law statement which shows that absolute zero cannot be reached 
since the rate of cooling vanishes as absolute zero is approached. 

B. The asymptotic relation between the internal and external temperature on the 

cold branch 

When the bath temperature tends to zero, the internal working fluid temperature has to 
follow. This becomes a linear relationship between T' and T c as T c tends to zero. 
Calculating the polarization at the end of the contact with the cold bath S2: 

o Q e q (S e 2 q ~ St q )x{l - y) - a 2 x{l/T b +y/r a ) 

(1 - xy) 

Assuming the relation Th = pT c as T c tends to zero, the exponents can be expanded to the 
first order to give: 

T c 1 - xy max + (uj a /uj b )px(y max - y) x(a 2 /r a )(y max - y) 

82 ~ £ 0^ + 1/2 i^y) (3 ' 8) 

Also, ^2 defines the internal temperature T' through the relation: S 2 = —\ tanh ( 2ksT' ) " 
Expanding the hyperbolic tangent, one gets: 

T , _ T 1 ~ xy max + p(uj a /uj b )x(y max - y) _ xu a (a 2 /r a ){y max - y) 

(1 - xy) (1 - xy) 

proving that T c and T' both tend asymptotically to zero. It should be noted that the term 
independent of T c depends on ui a , which also tends to zero as T c tends to zero ( (Eq. |3.6[ ) . 
Eq. ( |3.6| ) also shows that Q P Umum * r £ c j s a quadratic function of uo a) Cf. Fig. ^. 
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Eq. ( p.6|) represents an upper-bound to the rate of cooling. In order to determine how 
closely this limit be approached, a strategy of cooling must be devised, which re-optimizes 
the cooling power during the changing conditions when T c approaches zero. 



C. Optimal cooling strategy 

The goal is to follow an optimal cooling strategy, which exploits the properties of the 
equations and achieves the upper-bound for the rate of cooling, Qp. 
The properties of the equations employed are; 



• i: The derivative with respect to u a of Qp ( Eq. ( |3.4|) ), is extensive of order zero in 
the 'quartet' (u a , u b , T c , T h ). 

• ii: For the extensivity holds also for the 'doublets' {uo b , T h ) or {uo a , T c ). Scaling 
these variables by the same number, leaves Eq. ( |3.4j ) equal to zero, and the value of 

Q optimum j , i _ 

/ does not change. 

• iii: In spite of Qp being monotonic in u b} Q°^ timum i s independent of u b (and of T^), 
therefore Qp saturates as u) b is increased. 

From property (i) it follows, that once an optimal 'quartet '(o> a , u b , T c , T h ) is created, it 
is possible to cool optimally with a set of quartets, which are scaled by a decreasing set 
r n < 1, lim n ^oor n = 0. For this set the limit of the ratio ^ is a non zero constant. 

J- c 

Therefore in Eq. ( |3.5|) uo a and Tc are optimal leading to: 

^optimum -^(-^j Ij) ( ^ a,optimal \ T^cfiptimal /q in\ 

^ >^B 'optimal \ - 1 c, optimal J CUoIl ^ fc B p ( . ; 

In general, the hot bath temperature is constant, and the property (ii) is used to scale 
back the value of the optimal T^ to the bath temperature. As a result, the optimal high 
field is also scaled. 

Property (iii) will be exploited by changing only u b in the optimal quartets and checking 
for saturation. See Fig. [13] and the dashed curves of Fig. ||. Summarizing, for every 'quartet' 
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the upper-bound in Eq. ( |3.6|) can be reached. The details of the cooling strategy can be 
found in Appendix [E] 

1.5 i ■ 1 ■ 1 ■ 1 i ■ 1 ■ 1 ■ 1 
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FIG. 9. Left, (solid line): The optimal heat-fiow Q°P timum f or ^he heat pump 
as a function of T c . The fixed parameter values of the pair (a»j and T/J are, 
T h = 64.359, 51.49, 42.9082, 30.03555, 23.599, 17.153, 4.291, 0.8582, 0.1717, 0.01717, and accordingly 
u b = 150,120,100,70,55,40,10,2,0.4,0.04. The other constant parameter values are: a = 0.005, 
T c = 1, Th = 2. (dashed line): fixing ujj, = 3000 for every point. The other parameters are the 
same as the solid line. Right: The optimal heat-flows multiplied by the corresponding T c as a 
function of u> a . The optimal time is constant for the chosen parameters; Tpptimal = 

0.885 for the 

solid curves, and 0.825 for the dashed curves. 

Fig. [K] shows that the cooling strategy ( Tables [V] and |VT| ) can approach the upper 
bound leading to a linear relation of the optimal cooling power with temperature. With 
respect to the fields the optimal strategy leads to a decrease of the field uo a which is in 
contact with the cold bath. This causes the internal temperature of the TLS T' to be lower 
than the cold bath temperature T c . On the hot side the optimal solution requires as large 
an energy separation as possible u a — > oo but this effect saturates. 
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FIG. 10. Left: The optimal cooling rate as a function of the cold bath temperature T c , com- 
pared with the upper-bound for the cooling rate. Right: The entropy production during cooling 
shown for the case with friction (upper line, circles) and without friction (lower line, squares). The 
common parameters for all three figures are: r a = t& = 0.01, a = 0.005, T c = 1, = 2. 

The linear relation of the cooling rate with T c leads to a constant asymptotic entropy 
production as can be seen in the right of Fig. [TTj ( Cf . Appendix C) . 




IV. CONCLUSION 

The detailed study of the four stroke discrete heat engine with internal friction serves 
as a source of insight on the performance of refrigerators at temperatures which are very 
close to absolute zero. The next step is to find out if the behavior of the specific heat 
pump described in the study can be generalized. A comparison with other systems studied 
indicates that the conclusions drawn from the model are generic. As a heat engine the model 
shows the generic behavior of maximum power as a function of control parameters found 
in finite time thermodynamics |3p||. This is despite the fact that the heat transfer laws 
in the microscopic model of the working fluid are different from the macroscopic laws such 



as the Newtonian heat transfer law [16]. When operated as a heat pump with friction, the 



present model shows the universal behavior observed for commercial chillers ]27] caused by 
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a tradeoff between allocating resources to the 'adiabats' or to the heat exchange branches. 

Another question is whether the linear scaling of the optimal cooling power at low cold 
bath temperatures is a universal phenomenon. For low temperatures the results of the 
present model can be extended to a working fluid consisting of an ensemble of harmonic 
oscillators or any N-level systems. This is because at the limit of absolute zero only the two 
lowest energy levels are relevant. When examined, other models with different operating 
cycles show an identical behavior. For example the continuous model of a quantum heat 
engine [Tj| based on reversing the operation of a laser shows this linear scaling phenomena. 
Another example is the Ericsson refrigeration cycle Cf. Eq. (23) in the study of Chen et al 



25] which shows the same asymptotic linear relationship. 

A point of concern is the dependence of the heat transfer laws on temperature when 
absolute zero is approached. The kinetic parameters k± and fcf represent an individual 
coupling of the two-level-system to the bath. Considering coefficients derived from gas 
phase collisions they settle to a constant asymptotic value as the temperature is lowered 



|2T| . The reason is that the slow approach velocity is compensated by the increase in the 



thermal De-Broglie wavelength. 

There has been an ongoing interest in the meaning of the third law of thermodynamics 
|30H36|. The issue at stake has been: is the third law an independent postulate or it is a 
consequence of the second law and the vanishing of the heat capacity. This study presents 
a dynamical interpretation of the third law. The absolute temperature cannot be reached 
because the maximum rate of cooling vanishes linearly at least with temperature. 
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APPENDIX A: ANALYSIS FOR THE 'MOVING' CYCLES. 



Insight into the origin of the behavior of the 'moving' cycles is seen in Fig. 11, where the 
polarizations Si, S2 are shown as monotonically decreasing functions of the time allocation 
on the cold bath. However, the envelope of Si for maximal power, namely for maximal S1-S2 
is worth noticing. It is a decreasing function for short cycle times, achieves a minimum 
at ro, and starts to increase for r > To- Thus it is responsible for shifting the cycles to 
smaller polarization for short cycle times, and for the change of that trend for larger cycle 
times. The envelope of S2 for maximal S1-S2 is also a monotonically decreasing function of 
r c , or equivalently of r, supporting the increase of S1-S2. The figure also shows, that for a 
short time allocation both Si and S2 are close to the equilibrium polarization S^ 9 , When not 
enough time is allocated on the hot bath both the polarizations Si and S2 approach S^ g . 
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-0.48 



FIG. 11. Comparison between the polarizations Si and S2 as a function of r c , for ten different 
r values, 0.06, 0.08, 0.1, 0.25, 0.5, 0.75, 1 1.25, 1.5 and 2. The solid curves are Si while the dashed 
curves are S2. Superimposed are the values of Si and S2 for the maximal Si - S2. 
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APPENDIX B: THE COMPUTATION OF r 



The computation of tq Eq. ( 2.39 ) is not sufficient since it gives only the relation be- 
tween the times spent on the cold and hot branches for zero work. The natural addi- 
tional requirement is to seek for the optimal allocations, r Cj o and Th,o using Eq. ( 2.52 ): 

TO = Tcfl + T hfi + T a + T b 

Denoting by x and y the corresponding x and y values defined in Eq. ( f2.19| ), the 
following two equations for xo and yo are obtained: 

y° = Tr -r) - Rr (B1) 

and 

r c X ((1 — 2/ ) (1 - VO Xmax)) = I\ Vo ((1 ~ X ) (Xmax ~ X )) (B2) 

Where R is defined as: 



a 2 u a (l/r a + l/n) 
(uj b - uj a ) (Sl Q - Sc q + <J 2 /r a 



U UJ a K L/ , a T 1/ l b ) , . 

K ~ 7~. . . \ / Q e Q Qeq I -9 /- \ y^ 6 ) 



and x max was defined in Eq. ( |2.29[) as: 



(Sf -S?) + <T*/T a { ' 

The quadratic equation to be solved for x is, 

AA x\ + BB x + CC = (B5) 

Where AA = T h (1 + R) BB = - (T h ((1 + R) (x ma ,.-R) + x ma ,) + T c (1 + R - x max )) 
and CC = I\ (x max -R) x max 

APPENDIX C: ENTROPY PRODUCTION. 
(1) Heat Engine. 

A«S^ el = - (Qab/Th + Q CD /T C ) (CI) 
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Or from Table (|) 

A^ el = (ujTc-ut/T^-Sz) + °^{l/r a + l/r 6 ) (C2) 

The entropy production results are shown in Fig. |1^. The left figure shows AS U with 
increasing friction. The middle figure shows the corresponding cycles, while the right figure 
shows the corresponding power values. 




\ CO \ 



FIG. 12. Left: The entropy production of the heat engine as a function of the time spent on the 
cold branch for the fixed values of uj a = 1794, u>b = 4238, T c = 500, = 2500, r c = 1, Th = 2 and 
r = 0.1. Middle: The corresponding cycles. Right: The corresponding powers. Seven cases are 
shown. Case '00' is the frictionless case, when the times spent on the 'adiabats', are zero, case '01' 
is the frictionless case, when the times spent on the 'adiabats', r a and are different from zero and 
equal: 0.01. The other five cases are with increasing friction, when also r a = t\, = 0.01, whereas 
the different friction coefficients a are: for plot 1: a = 0.003, for plot 2: a = 0.004, for plot 3: 
a = 0.005, for plot 4: a = 0.006 and for plot 5: a = 0.007. 

The reciprocal behavior of the entropy production and the power is clear from Fig. [T^. 
One also observes, that for the given cycle time the 'free' time for the cycles with increasing 
a becomes more restricted. This follows from the dependence of T c ^ min on o. See also Fig.( 
1) 

Introducing Eq. (|2.28|) into Eq. (|C2|) . The entropy production becomes, 

AS: ylel = {u a /T c -u b /T h ){S<?-S?)F{x,y) + (C3) 

where 
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A<S"i = ^J^h^ (fr(l-^)(l/ra + yM) + ^(l-y)(x/r a + l/r b ) ) (C4) 



Notice, that A«S" 1 is always positive. For a = Eq. ( |C4| ) reduces to the frictionless results 

zg. 

(2) Heat Pump 

The entropy production for the heat pump becomes: 



2 771/ n/^ 1 z 1 a\ w a 1 A y A . . 

aF{x,y)\—- (— + -) + — (- + — ) > C5 

[T h l- x r a r h T c 1 - y n r a J 

The asymptotic entropy production as T c tends to zero can be calculated leading to 
A<S? e/ = F(x,y) \{u} b /{puj a )){l- P {uj a /uj b )) 2 + 



a 



2 



^(1-*)^ + ^ + t^y) (1/Tt + vM ) ] (C6) 

Since T h = pT c , the r.h.s. of Eq. (|C6| ) tends to a constant, for each term depends on 
the constant ratios (o^/T^), (u a /T c ) or on their ratio. This result is demonstrated on the 



right side of Fig. 10 



The optimization with respect to time allocation has the same result as for the heat 
engine. Therefore, only optimization with respect to the fields are presented; 

Equating to zero the derivatives with respect to x an y of the entropy production, one 
gets two similar equation to the total work derivatives: 

( I- - V X max ) + ^ / UJa_\ + 1-3/ > (C7) 



(u a /T c - uj b T h ) \2k B Tj (4:k B T ( 



(AS*« + a 2 /r a ) cosh 2 (-^-) + X >0 (C8) 

{u) a /T c - u b T h ) \2k B T h J (4 k B T h ) 

Where AS eq , is S e h q - S e c q . 

Eqs. flU7|) and ( |C8| ) show that the entropy production is a monotonic function in the 
allowed range, namely, for 
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To conclude the entropy production has a minimum value: AS^ in , will be 

A5^ n = |V(1/t + l/n) (CIO) 



c 



obtained on the boundary of the range. 

APPENDIX D: THE TOTAL WORK DONE ON THE SYSTEM FOR THE HEAT 

PUMP 

The total work done on the system becomes, 

W^ e3 = (oo b - 0J a )(S 2 - Si) + a 2 uo b (l/r a + l/n) (Dl) 

or 

W° cy n le3 = (u b - u a ){S e 2 q - St q )F{x, y) + W a3 (D2) 

where 

2 

W ct3 = 7-^ r(o; 6 (l-y)(l/r a + x/r b ) + u a {l - x){y / r a + l/r 6 )) 

(1 - xy) 

= ^(x, w ){J5-(I + 1) + -^(1 + »)} (D3) 

Eq. ( pi| ) can be interpreted as the work done on the working fluid see (Cf. Fig. 0), as 
the sum of three positive areas, (u b — co , a )(S , 2 — Si), a 2 u b (l/r a ) and a 2 u b (l/T b ) with the 
corresponding corners, D^B^A 1 , B,B 1 ,S 2 ,S 3 and A 1 ,A,S 4 ,Si. 

APPENDIX E: THE OPTIMAL COOLING STRATEGY CLOSE TO THE 
ABSOLUTE ZERO TEMPERATURE 

The first step in the cooling strategy is to create the first optimal quartet; 

• (0) The systems external parameters a, r a , r b , T c and 1^ are set. 
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(1) A decreasing set of Ub is chosen. 

(2) A constant ratio (p) for Th/T c , is chosen which is the ratio of the initial bath 
temperatures. 

(3) For the above chosen values, the optimal values of ou a , T c , t, and its optimal 
allocations between the branches to give maximal are found for each u>b in the set 
in (1) , by solving numerically the following additional equation to Eq. (|3.4j ), with the 
condition that T h = p-T c : 

dQ F F(x,y) f uj a \ 2 I 1 u b 



dT c At k B V T C J \ cosh 2 "° T pu a cosh 2 



(El) 



2 k B T c 2p k B T c 



The above strategy causes the decrease of Th together with the T c . Nevertheless according 
to (ii) above, the doublet u>b and Th can be rescaled to increase Th back to its original value. 
The solid curves of Fig. |9| are optimal in the in the above described sense. Increasing only 
the value of Uf, in the optimal quartet according to point (iii), leads to larger values of the 
cooling rate, but eventually the increase of Qp will slow down and saturate. See Fig. O 
and the dashed curves of Fig. |S|. 

Fig. ^ represents the saturation phenomenon on Ub- Three points from Fig. ^ are 
chosen, and all parameters are fixed, except Ub, which is allowed to increase. 
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FIG. 13. Left: The optimal heat-flow for the heat pump as a function of Uf,, showing 
the saturation phenomenon. The fixed parameter values are: for triangles; Th = 64.5725, 
T c = 12.9145, uj a = 11.9233, starting with oo b = 150, for squares; T h = 42.90815, 
T c = 8.58168, oo a = 7.94986, starting with u b = 100, for circles; T h = 23.59905, T c = 4.71981, 
Lo a = 4.37247, starting with cob = 55. The common parameters for all three figures are: 
t = Tb = 0.01, a = 0.005, T c = 1, Th = 2. Right: The optimal heat-flows as a function of 
T maxi the time at which the optimum is achieved. The fixed parameter values are the same, as on 
the left. We note, that the optimal time is becoming constant only at saturation. 



In order to approach the upper-bound for Qp in Eq. (|3.6|), a decreasing set of u a / T c is 
created, achieved in an optimal way: 

First step: After having an optimal 'quartet', T c and T^, are fixed. Then, by lowering 
u b , one finds the corresponding optimal uj a values. This procedure is checked globally, by 
also iterating the time allocations. The results of a typical example are shown in Table [V]. 

Second step: Using again the property of extensivity, the cooling will be achieved by 
multiplying the rows of Table by a decreasing sequence, e.g. by 2 _n for the n-th row. 
Table |VI] describes the cooling strategy, checking also the non-divergence of the entropy 
production both for the frictionless case and the case with friction. The results are also 



summarized in Fig. 10 
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Table [V| demonstrates, that the procedure shifts down to the Carnot bound. The ratio 
R = jp-/ ^a. was computed showing only small changes. 

TABLE V. First step. Starting from an optimal quartet, the procedure creates for a given 
decreasing set of cj&-s, a decreasing set of of u> a -s for fixed bath temperatures. 



T c 


T h 




, .optimal 


^optimal 




R 


Qf 


0.0025 


50 


60 


1.370(-3) 


0.5392 


1.2 


2.226 


5.812(-5) 


0.0025 


50 


55 


1.273(-3) 


0.5090 


1.1 


2.161 


5.013(-5) 


0.0025 


50 


50 


1.164(-3) 


0.4653 


1 


2.149 


4.241(-5) 


0.0025 


50 


45 


1.051(-3) 


0.4205 


0.9 


2.140 


3.50549(-5) 


0.0025 


50 


40 


9.320(-4) 


0.3728 


0.8 


2.146 


2.826(-5) 


0.0025 


50 


35 


8.250(-4) 


0.3300 


0.7 


2.121 


2.182(-5) 


0.0025 


50 


30 


6.985(-4) 


0.2794 


0.6 


2.147 


1.613(-5) 



TABLE VI. A procedure to get an optimal set of pairs of u a , T c where their ratio tends to 
zero. Th = 50 for every cold bath temperature, T c . The index 'fl' stands for the frictionless case, 
and Q U p denotes the upper-bound for Qp. 



T c 


0J b 


, .optimal 
u a 


AS U 


AS uJl 


Qf 


Of 


0.0025 


60 


1.370(-3) 


0.0333353 


0.0285 


5.812(-5) 


6.084(-5) 


0.00125 


55 


6.365(-4) 


0.0285075 


0.02406 


2.457(-5) 


2.626(-5) 


0.000625 


50 


2.91(-4) 


0.0242662 


0.02021 


1.039(-5) 


1.098(-5) 


0.0003125 


45 


1.3138(-4) 


0.020338 


0.01669 


4.2939(-6) 


4.467(-6) 


0.00015625 


40 


5.825(-5) 


0.016788 


0.01354 


1.7239(-6) 


1.759(-6) 


0.0000781 


35 


2.578(-5) 


0.013210 


0.010357 


6.6772(-7) 


6.888(-7) 


0.0000391 


30 


1.0914(-5) 


0.010356 


0.007915 


2.4673(-7) 


2.468(-7) 
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